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Abstract 



Within the composite operator formalism we derive a formula for the pion de- 
cay constant f w , as defined directly from the residue at the pion pole of the meson 
propagator, rather than from the matrix element of the axial current. The calcula- 
tion is performed under some simplifying assumptions, and we verify the complete 
consistency with soft-pion results, in particular with the Adler-Dashen relation. The 
formula one obtains for (the pole-defined) f% differs from the previous Pagels-Stokar 
expression by an additive term, and it still provides f% in terms of the quark self- 
energy. We make some numerical estimates leading to (30 -j- 40)% deviation for f% 
with respect to the Pagels-Stokar formula. 
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1. Introduction 

In quantum chromodynamics the pion, as well-known, is a pseudo Nambu-Goldstone 
boson associated with the spontaneous breaking of chiral symmetry. Its decay constant f n 
plays a key dynamical role in the chiral symmetry breaking mechanism of QCD, and anal- 
ogous quantities appear in other theories which use similar paradigms, such as electroweak 
symmetry breaking through a new strong sector. In 1979 Pagels and Stokar proposed an 
approximate expression for the calculation of f„ Jl|. Their derivation used a sum rule due 
to Jackiw and Johnson 0, plus additional assumptions within the so-called dynamical 
perturbation theory, and allowed for an approximate expression for f n , as defined from 
the matrix element of the axial current, in terms of the quark self-energy. The formula is 
of great utility in QCD and also in theories derived from the old technicolour concept. 

There has been a vast literature on the Pagels-Stokar (PS) formula with the general con- 
clusion that it leads to a sensible result for f n , within uncertainties not easily controllable 
in view of the theoretical approximations present in the derivation ]3[]. The Pagels-Stokar 
expression for the pion decay constant is 



where d(r) is the dimension of the quark colour representation (d(r) = 3 in QCD) and E 
is the dynamical quark self-energy in the chiral limit. 

We shall present below a new formula for which shares with the PS formula the 
advantage of only depending on the self-energy S , and that we derive within the com- 
posite operator formalism developed in ref. |4| and as modified in ref. ||. Within such 
schemes one could think of two different calculations of the pion decay constant. One can 
evaluate the coupling of the pseudo-Goldstone to the axial- vector current. Alternatively, 
one can directly evaluate the residue at the pion pole of the meson propagator. These 
procedures correspond to different definitions, but only the second one agrees with the 
current-algebra and soft-pion results || (in particular the Adler-Dashen relation). We 
shall follow the second procedure, that is calculating the residue at the pion pole. Within 
the composite operator formalism, for a first approximate understanding, we discuss here 
the so-called "rigid case", in which the presumably small logarithmic corrections com- 
ing from the renormalization-group analysis are neglected. For the asymptotic behaviour 
of the self-energy function we choose the one dictated by the operator product expansion 
(OPE) [Q] up to a logarithm which we neglect. The motivation for such approximate study 




, 2 ^S (fc 2 ) 



(1) 
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is mainly that in this way it is possible to derive analitically the complete expression for 
the effective action at two fermion loops. Furthermore a phenomenological analysis we 
have previously performed of the pseudoscalar masses || in this contest did not show any 
major inadequacy of such a treatment. Our expression for f 2 is 
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The first term in Eq. (2) is fp S of Eq. (1). By writing f 2 = f PS (1 + 5 2 ), to get an 
evaluation of 5 we can go back to two different Ansatz for S we had used in the past 
to study low energy QCD || |7j. We find for 5 2 values such as 0.35 and 0.37, suggesting 
that the correction to the Pagels-Stokar expression is presumably large (30 4- 40)% and 
presumably rather insensitive to the modification of the self-energy shape, provided the 
ultra-violet behaviour in k 2 is roughly maintained. 



2. The effective action 



We start from the effective Euclidean action for the composite operator formalism 



r(E) 



-Tr In 



ST2 
5S 



+ Tr 



ST2 
5S 



r 2 (S) + counterterms 



(3) 



where S 1 



(ip — m), m is the bare quark mass matrix, ^(S) is the sum of all two- 
particle irreducible vacuum diagrams with fermionic propagator S and S = — Eq. 
(3) is the modification of the effective action of Cornwall, Jackiw and Tomboulis 0] which 
was introduced in ref. || to account for the correct stability properties of the theory. At 
two- loop level T 2 = ^Tr(SAS), where A is the gauge boson propagator, so that S = — AS, 
Tr [£T2/5S] = 2T2, and one can rewrite Eq. (3) in terms of S 



— Tr In 
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(4) 



Here the variable S plays the role of a dynamical variable. At the minimum of the 
functional action, that is when the Schwinger-Dyson equation is satisfied, S is nothing 
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but the fermion self-energy . A parametrization for S, employed in H, was 

£ = ( s + z 75 p)/(A;) = S S + 2 75 S P (5) 

with a suitable Ansatz for f(k), and with s and p scalar and pseudoscalar constant fields 
respectively. 

3. The effective potential 

The effective potential one obtains from Eq. (4) (see ref. ||) is 

T 8vr 2 d(r) f d A k r _ „ „ 

V=— = ; K ~' / t — rrtr S s n fe S s + 53„n fc SJ - 

^ 3C 2 (r)^ 2 i (2vr) 4 1 s fe s p * pJ 

-d(r) Tr In [ijfe - (m + S s ) - r^EJ + 5Z tr(ms) (6) 



where C^r) is the quadratic Casimir of the fermion colour representation (for SU(3) C 
C 2 = 4/3) and S s = \ a s a f(k)/V2, S p = \ aPa f(k)/V2, m = \ a m a /y/2 (a = 0, • • - , 8, 
Ao = y2/3, Aj = Gell-Mann matrices, z = 1, • • • , 8). Furthermore has a divergent piece 
to compensate the leading divergence proportional to tr(ms) in the logarithmic term. For 
a quark of mass m the effective potential is 



V(s,p,m) = -d(r)cj—^[E s n k E s + E p D k E p 



d A k 

-2d(r) [ 44 ^ \k 2 + (m + S s ) 2 + £ 2 1 + 5Zms (7) 

where we have defined c = 2n 2 /g 2 . In ref. 0] after fixing so as to cancel the lead- 
ing divergence proportional to ms in the logarithm, we had imposed the normalization 
condition 

1 (8) 



1 dV 
hm — 



m^o m d{ipip) 
or, with (ijnf)) = (d(r)M 3 /2rr 2 ) cs 

1 dV 



extr 



lim 

m-»o as 



extr 



rf(r)M3 c 



where M is a momentum scale for the self energy and s is the extremum of the effective 
potential. The extrema of the effective potential in the massless case m = depend only 
on c. Therefore Eq. (9) becomes, in this case, an equation for c and M is left undetermined. 
This is nothing but the usual dimensional transmutation. The numerical values for c and 
So, the minimum of the effective potential in the massless case, are obtained once one has 
fixed the Ansatz for f(k). 
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With Sq(A;) = s f(k) we shall write in general 



5Z = d(r) 



dV 



d 4 k 
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The gap equation, from : -— = 0, is 

ds 



d{r) 
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d A k 
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d 4 A; (m + E s )E s 



(2tt) 4 A; 2 + (m + E s ) s 



+ m5Z = 



(11) 



where E s = sf{k) and s is the value at the minimum in the presence of the bare mass. 

Let us now turn to the effective action. The fields s and p depend in this case on the 
space coordinates and we shall use the Weyl symmetrization prescription 



E = (s + n 5 p)f(k) -> - [s(x) + * 7 5P(x), /(k)]_ 



(12) 



We are interested in oscillations around the minimum of the effective potential, so we 
introduce 



X(x) = s(x) - s, tt(x) = p(x) - p = p(x) 
v(x) = x(x) + i7 5 7r(x) 
S(Jfe) = ik- ((m + S s (A;)) 

The Tr In term in Eq. (4), which we denote as Ti og , becomes 



T log = -Trln 



ik — m — E 



= -Tr In 



s- 1 --^),/^)]. 



(13) 



(14) 



As we are interested in the 2-points function, we expand to second order in v(x) and after 
some calculation we obtain for the Fourier transform of Ti og in Eq. (14) 

„ j _ r d 4 k , r, 9 / ^ /1N \2l f d 4 k 
Ti og = d(r){-2nj — — trln^ 2 +(m + E s (A;) 
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where 



1 



V^fca) = -[/(*!) + /(fc 2 )] 



(16) 
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For T 2 we obtain (working in Landau's gauge) 

d 4 k 



T 2 = -d(r) ctt J 



(2vr) 4 
d 4 k 



tr 



S s D fe S s 



+ 1 |H - d{L)c< - q) {I Wr mUkf{k) ) <q) 



/(*)□*/(*) J tr(i X (0)) 
d 4 k 



+ (pseudoscalar <-> sca/ar) j 
For the counterterm one has 



(17) 



T ct = tr(ms)fi<5Z + tr(mx(0))6Z 



4. The improved expression for 



We note that each term in the effective action consists of a constant, a linear and a 
quadratic term in the fields. The constant term gives back the original potential Eq. (6) 
at the minimum. Such a term controls the normalization. The linear term vanishes by 
virtue of the gap equation, Eq. (11). The quadratic term stands up for the effective action 
up to the second order in the fields. In space-time coordinates 



Jd*xd*yJ-0- 4 e- i( l( x -y)n a (x) 



V2 
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d 4 k 
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(19) 



From 



5 2 T 



57T a (x)5^(y) 
one finds for the Fourier transform of G~^(x — y) 



(20) 
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(21) 
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By using the gap equation to eliminate c we get, for a quark of mass m 

d 4 k 



d(r 
4 



rtr 



(27T) 4 

Ad{r) r d 4 k 
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(2tt) 4 k 2 + (m + S s ) 2 s"~ ^ 

We can eliminate 5Z by using the normalization condition, see Eq. (10), and the 
relation {ipip) = (d(r)M 3 /2ir 2 ) cs to obtain 

4m r d 4 k E (k) 



G-\q 2 ) = -- 2 m + d(r)\-— f 
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(23) 



Note that the second and third terms in Eq. (23) regularize each other in the ultraviolet 
In the limit of small quark masses, expanding in we find 

g-V) = (4-V-(? 2 -^<^> 



.V2s 0/ 



(24) 
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(25) 



In Minkowski metrics the propagator (24) has a pole at q 2 = m 2 = -2m(^) /F 2 , with 
residue (a/2sq/-F) 2 , where (ip^o = (d(r)M 3 /2h 2 )csq. The Adler-Dashen relation (which 



follows from the symmetries and current algebra) requires the identification f 2 



J , so 



that q 2 



ml 



-2m{ipip)o/ f 2 - The rescaling factor b relating the canonical field ip^ (with 



unit residue at the pole) to the field n, (p n = bn, is then b = —f w /y/2so, as indeed follows 
from current algebra and soft pions theorem (see ref. ||). Comparison of (25) with the 
Pagels-Stokar formula (1) leads to our new formula (2) of the introduction. 

To get a numerical insight into the problem we use the dynamical calculations in 
ref. [||, 0, where overall fits to low energy QCD were made on the basis of two alter- 
native Ansatz for £(&) = sf(k): a smooth Ansatz f(k) = M 3 /(M 2 + k 2 ) for which 



8 



the relevant parameters took values c = 0.554, s = —4.06, and a step-function Ansatz 
f(k) = M[6(M 2 - k 2 ) + (M 2 /k 2 )9(k 2 - M 2 )] for which one found c = 0.32, s = -2.69. 
Our new expression for has f 2 = fpsO- + $ 2 )-, where 8 follows from Eq. (25). We find 
5 2 = 0.347 in the case of the smooth Ansatz, and 8 2 = 0.376 for the step-function Ansatz. 
We want to remark that in the massless limit we consider, the correction 8 depends only on 
So (or c) and the shape of f(k), but not on M. In particular, because of the fact that the 
relevant contribution to the chiral symmetry breaking phenomenon comes from relatively 
short- distance effects, the corrections will depend mainly on the ultraviolet behaviour of 
the self-energy. That is, the correction does not depend on the fit we have made to low 
energy QCD. 

It therefore seems that: (i) the corrections are relevant with respect to the old Pagels- 
Stokar formula; (ii) the corrections do not seem to vary in a sensible way when varying 
the Ansatz for the self energy, at least within the Ansatz we have used. 

It is obvious that the next step one should take is to see whether the neglected cor- 
rections, which we know must be there, can modify the results. However, due to previous 
experience from the study of the pseudoscalar masses |J, we would not expect substantial 
changes in the overall picture of dynamical symmetry breaking. 

Finally we may stress that the new formula (2) we have obtained for /„., within the 
approximations made, does not require additional inputs beyond those already present in 
the Pagels-Stokar formula. 



This work has been carried out within the program Human Capital and Mobility 
(BBW/OFES 95.0200; CHRXCT 94-0579) 
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